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Abstract 

Equilibrium states in galactic dynamics can be described as station- 
ary solutions of the Vlasov-Poisson system, which is the non-relativistic 
case, or of the Vlasov-Einstein system, which is the relativistic case. 
To obtain spherically symmetric stationary solutions the distribution 
function of the particles (stars) on phase space is taken to be a func- 
tion <&(E,L) of the particle energy and angular momentum. We give 
a new condition on $ which guarantees that the resulting steady state 
has finite mass and compact support both for the non-relativistic and 
the relativistic case. The condition is local in the sense that only the 
asymptotic behaviour of <f> for E^E needs to be prescribed, where 
Eq is a cut-off energy above which no particles exist. 



1 Introduction 

Let / = f(t,x,v), telR, x,uGlR 3 be the density function on phase space 
of the stars in a galaxy. We assume that collisions among the stars are 
sufficiently rare to be neglected and that the stars interact only by the grav- 
itational field which they create collectively. In a non-relativistic situation 
/ obeys the Vlasov-Poisson system, in a general relativistic situation one 
obtains the Vlasov-Einstein system. We are interested in spherically sym- 
metric equilibrium solutions of these systems, i. e., / is independent of time 
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t, and f(Ax,Av) = f(x,v), t4gSO(3), x,-u£lR 3 . The Vlasov-Poisson system 
then takes the form 

v-V x f-U'--V v f = 0, (1.1) 

r 



4tt 



r 



U>( r )= s 2 p(s)ds, (1.2) 



p(r)=p(x) = J f(x,v)dv. (1-3) 

Here r = |x|, ' denotes derivative with respect to r, [7 denotes the gravita- 
tional potential of the system and p the spatial mass density induced by /; 
we assume that all the stars have mass one and note that due to spherical 
symmetry U and p depend only on r. The Vlasov-Einstein system takes the 
form 

:-V x /-Vl + «V--V„/ = 0, (1.4) 



VI + v' 

e~ 2A (2rA'-l) + l = 8vrr 2 p, (1.5) 
e~ 2A (2r/i' + l) -l = 8vrr 2 p, (1.6) 



p(r) = p(x) = J \Jl + v 2 f(x,v)dv, (1-7) 

p(r)=p(x)= J f(x,v)^L=. (1.8) 

Here v 2 = v-v, and • denotes the Euclidean scalar product on IR 3 . If x = 
(rsin6cos(j),rsm8sm(fi,rcos6) then the space time metric is given by 

ds 2 = -e 2lM dt 2 + e 2X dr 2 + r 2 ( d6 2 + sin 2 



p denotes the mass density and p the radial pressure. As to the choice 
of coordinates on phase space, which leads to the above form of the static, 
spherically symmetric Vlasov-Einstein system, we refer to [14]. As boundary 
conditions we require asymptotic flatness, i. e., 

lim p(r)= lim A(r) = 0, (1-9) 

r^oo r^oo 

and a regular center, i. e., 

A(0) = 0. (1.10) 
For the Vlasov-Poisson system we require 

lim[/(r) = 0, (1.11) 
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and we will also need the radial pressure p in the non-relativistic case, which 
is defined as 

p(r)=p(x)=J(j^\ f(x,v)dv, (1.12) 

There are essentially two approaches to construct solutions of these systems. 
The first is to observe that there exist invariants of the particle motion, 
namely the particle energy E = E(x,v) and angular momentum squared L = 
L(x,v); for the definition of these quantities cf. Q2.1| ) and ( p^ ) respectively. 
The ansatz 

f(x,v) = $(E,L) (1.13) 

with some prescribed function <3? automatically satisfies the Vlasov equation, 
and it remains to solve the field equation(s) with the ansatz for / substituted 
into the definitions for p and p; these quantities become functionals of U or 
jjL since E depends on U or fj, respectively. The main problem then is to 
show that the resulting steady state has finite (ADM) mass and compact 
support. In Q this was done in the non-relativistic case for the so-called 
polytropic ansatz 

$(E,L) = (E -E) k + L l , 

where Eq is some constant, (•)+ denotes the positive part, and k> —1, l> 
— 1, k+l + 1/2 > 0, k <3l + 7/2. In [10, 16] an existence result was estab- 



lished for the Vlasov-Einstein system, exploiting the fact that the Vlasov- 
Poisson system is the limit of the Vlasov-Einstein system as the speed of 
light tends to infinity, cf. p5|| . (In the present paper the speed of light is 
set to unity.) This perturbation argument does not give good control over 
the class of models obtained, a fact which motivates the search for a better 
method. It should be noted that in order to obtain a steady state with finite 
mass <3? must vanish for energy values larger than some cut-off energy Eq. 

The second approach is to define an energy- Casimir functional which has 
the property that its critical points are steady states, and then show that 
this functional has a minimizer over a certain set of phase space densities 
/. This approach was used in Q for the Vlasov- Poisson system. It has the 
advantage that it provides a certain nonlinear dynamical stability property 
of the steady state obtained, and the resulting steady states are more general 
than the polytropic ones: Only certain growth conditions on $ need to 
be prescribed. Nevertheless, the latter approach also makes use of global 
properties of the function i. e., of growth conditions for E close to Eq 
and close to — oo. 

The present paper follows the first approach, but we give a new charac- 
terization of <&'s which lead to finite mass and compact support both in the 
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non-relativistic and in the relativistic case: Except for some mild technical 
assumptions we only require that 

$(E,L) = c(E - E) k + L l + O({E - E) k + +s )L l as E^E 

where 

11 3 

k>-l, l>~^, k + l + ->0, k<l + -. 

Thus the characterization is purely local: Only the asymptotic behaviour at 
E = Eq needs to be controlled. Since such steady states together with their 
features of finite mass and compact support persist under perturbations of 
the ansatz function <3? as long as the form of the asymptotic expansion at 
Eq is preserved, they might be called structurally stable. Another impor- 
tant point is that we obtain finite ADM mass and compact support in the 
relativistic case directly, i. e., without a perturbation argument as was used 
in |l(], [Hf ]. The assumptions are much more transparent than the smallness 
conditions required when starting from the Newtonian limit. 

The main idea of the argument is as follows: If i?€]0,oo] denotes the 
radius of the support of a steady state then one can show that Eq — U(r) — > 
or Eq — e^( r ) — > as r —* R-. One needs to show that R<oo, and this is done 
by expanding all the relevant quantities is terms of Eo — U(r) or Eo — e^ r \ 
obtaining detailed estimates for the behaviour of the solution of the field 
equations as R— . This approach was strongly motivated by the paper 
of Makino[Q where a technique of this kind is used to prove finite mass and 
compact support for relativistic, spherically symmetric stellar models, i. e., 
for steady states of the Euler-Einstein system. The connection between 
the two situations is as follows. If a spherically symmetric steady state 
solution of the Vlasov-Einstein system corresponds to a choice of which 
only depends on E, then the energy density p and the pressure p define a 
solution of the Euler-Einstein system describing a self-gravitating perfect 
fluid. (A similar relation between kinetic and fluid models holds in the non- 
relativistic case.) The asymptotic behaviour of <!> which is important in our 
results corresponds to the way in which p depends on p (equation of state) 
in the limit p^> 0. 

Our paper proceeds as follows: In the next section we derive the reduced 
problems which are obtained by substituting the ansatz for / into p and p 
in the field equation(s). We also show that a cut-off energy Eq is necessary 
in order to obtain finite mass and compact support. This shows that ex- 
cept for the form of the dependence on L our ansatz is quite general. The 
main result is then stated and proven in Section 3. In a last section we 
consider steady states which appear in the astrophysics literature and show 
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that our result applies to most of these and proves that these steady states 
have the physically very desirable properties of finite mass and compact sup- 
port. Clearly, there are some polytropic steady states with finite mass and 
compact support which are not covered by our result: For the polytropes 
one needs /c<3/ + 7/2 whereas we require k <l + 3/2. In the last section 
we also comment on this discrepancy from the viewpoint of the relation to 
static self-gravitating fluid bodies and from the viewpoint of recent stability 
results for the Vlasov-Poisson system; cf. ^] . 

To conclude this introduction we mention some further references which 
are related to the present paper. Steady states of the Vlasov-Poisson system 
with axial symmetry are constructed in [11|. Spherically symmetric steady 
states with a vacuum region at the center are constructed in [12] both for 
the Vlasov-Poisson and the Vlasov-Einstein system. Among these there are 
examples of relativistic steady states violating Jeans' Theorem, which holds 
for the Vlasov-Poisson system and says that all spherically symmetric steady 



states are obtained by an ansatz of the form ( 1.13 ), cf. (2Cj. In |T3| steady 
states of the Vlasov-Poisson system are constructed where the matter is 
concentrated on a plane, using the variational approach mentioned earlier. 
Concerning the initial value problem for the time dependent systems we 
mention @, ^, 19 1, where global existence of classical solutions to the Vlasov- 
Poisson system is established for general data. For the Vlasov-Einstein 
system the existence theory for the initial value problem is far less complete. 
We mention [14, |17| for the spherically symmetric, asymptotically flat case 
which is of interest here. 



2 The reduced problem 

Consider first the Vlasov-Poisson system. The particle energy and angular 
momentum squared 

E = E(x,v) = ^v 2 + U(r), L = L(x,v) = \xxv\ 2 (2.1) 

are conserved along particle trajectories, i. e., are constant along solutions 
of the characteristic system 

x 

x = v, v = —U'(r) — 
r 

of (1.1); r = \x\. If we make the ansatz ( 1.13 ) then upon substituting ( |1.13| ) 
and ( |2.1[ ) into Q1.3| ), p becomes a functional of U, and it remains to solve 
the resulting nonlinear Poisson equation; for the moment we only require $ 
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to be a nonnegative, measurable function. Similarly, for the Vlasov-Einstein 
system one finds that 



E = E{x,v) = e^ r) Vl + v 2 , L = L(x,v) = \xxv\ 2 (2.2) 
are constant along solutions of 



p x = — yl + v 2 p!(r)- 



1.6) where p and p become 
into (O) and C|). It 



and it remains to solve the system ( |1.5| ) and ( 
functionals of /i upon substituting fll.l3| ) and ( 2.2 
is a simple computation to see that in the non-relativistic case substituting 
(pit ) into (fT|) yields 

. . 2TT f°° [2r 2 (E-U(r)) dLdE 

p(r) = — / $(E,L)—== — — . (2.3 

HK) r 2 Ju(r)Jo V W2(E-U(r)-L/2r 2 ) V ; 

In the relativistic case we obtain 



2vr , r {l£l < E-dLdE 



p{r) = ^ e -Mr) I , <S>(E,L)- 

r J&VJo ^E 2 -e 2 ^ r ){l + L/r 2 ) 

(2.4) 

and 

2tt , s r°° /-rWe-^M-i) i 

p ( r )= e -Mr) / / ${E,L)\ E 2 -e 2 ^)(l + L/r 2 )dLdE. 

H J e K-r) Jo V 

(2-5) 

Before proceeding further we demonstrate that <3? has to vanish for large 
values of E if the resulting steady state is to have finite mass, i. e., if the 
quantity 

M= I p(x)dx (2.6) 



is finite; M is the total mass or the total ADM mass of the steady state 
respectively. Here and below it will be useful to recall that the solution of 
( |T~5|) satisfying ( |1.10|) is given by 

n,, 2m(r) , 

e~ 2X = l ^ (2.7) 

r 

where 

rr 

m (r)=4vr / s 2 p{s)ds, (2.8) 
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at least as long as 2m(r) <r. Using ( |2.7| ) one can rewrite ( |1.6| ) in the form 

„2A 



/U '(r)=e 2A ( —^p-+A-Krp{r) ) . (2.9) 



Theorem 2.1 Lei 3> : IR 2 — > [0,oo[ be measurable. 

(a) Let (f,U) be a solution of ( fill) , ftlll) , fiDl) in the sense that f(x,v) = 
&(E,L) and U G C 1 (]0,oo[) solves ^TJfy with substituted in. Let 
M<oo. Then Uoo = lim r ^ 00 ?7(r) < oo, and &(E,L)=0 a. e. for E> 
Uco, L>0. 

(b) Let (f,X,p) be a solution of flLh)), fcTh] ), flT7\), ftTl[) i n the 
sense that f{x,v) = <fr(E,L) and A, p G C^QO, oof) solve (TJ\), with 
($■40 , $-°V substituted in. LetM<oo. Then fioo = lim r _ )OC ^(r) < oo, 
and <5>(E,L) = a. e. for E>e fl °°, L>0. 

Proof : Let us first consider the non-relativistic case. Since 

0<U'(r)<-^, r>0, 
U is increasing and has a finite limit as r— >oo. By (|2.3| ), 

,oc roc f 2r 2 (E-U(r)) dLdEdr 

M = 8vr 2 / / / $(E,L)- 



Ju(r)J0 ' ' y/2(E-U(r)-L/2r' 2 ) 

> 8^ 2 / / $(£,£) / - . == dLdE. 

JUaoJo J^L/(2(E-U 00 )) y/2{E-U{r)) 

The integral with respect to r in the latter expression is infinite for any 
E> Uoo and L>0, which implies that <3? has to vanish for such arguments. 

Let us now consider the relativistic case. Since by ( |2.S| ) fi is increas- 
ing, the limit /Xqo = lim r _ +0O //(r) G] — oo,oo] exists. Since m(r) < M and thus 
g 2A(r) < 2 for r > 4M, and p< p, we conclude that 

M'W < 2 + 47rr P( r )^ > r ^ 4M ' 

and 

u(r)< J u(4M)+2 / -^-dr + 4vr / s 2 p(s)ds < oo, r > 4M 

,/4M H 4M J 4 M 
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which proves that /ioo <oo. Using (2.4) we obtain 

poo poo pr 2 (E 2 e~ 2 f J -( r ^-l) T? 2 rlTJT?rl r 

M = 8n 2 e- 2 ^) / / $(E,L)- L ' lLt1L(h 



ertr)Jo ' ^jE 2 -e 2 ^(l + L/r 2 ) 
>8ir 2 e- 2 ^ / / E 2 $(E,L) / dTAP, 

JeMoo Jo J^L/(E 2 e- 2 ^-l) \/ E 2 - e 2 ^ r ) 

and the assertion follows as in the non-relativistic case. □ 

Although this will become important only in the next section where we 
prove our main result we restrict the ansatz in the following way: 

f(x,v) = (j ) (E)L l (2.10) 

with l>—l/2 and <j) measurable. Observe now that in the non-relativistic 
case E can attain any real value, whereas in the relativistic case E > 0. In 
order not to introduce unnecessary assumptions this needs to be reflected 
in the choice of domain of </>, which we denote by Ji^njOof with 



E n 



-oo in the non-relativistic case, 
in the relativistic case. 



Since we are interested only in steady states with finite mass we assume 
that (p(E) = for all energies E larger than some given Eq > E m i n . It will be 
necessary to have some information on the functional dependence of p and 
p on U or [i respectively: 

Lemma 2.2 Let <p :]E m [ n ,oo[-^>TR be measurable, Eo>E m [ n , and k>—\ 
such that on every compact subset K c]E m [ n ,oo[ there exists C>0 such that 



Define 



0<<p(E)<C(E -E) k + , EeK. 

POO 

g m (u) = cj)(E)(E-u) m dE, ue]£ min ,oo[, 

J u 

POO 

h m (u)= / 4>(E)(E 2 — u 2 ) m dE, ue].E m i n ,oo[. 

If m> — 1 and k + m + 1 > then g m ,h m £ C(]E m i n , oo[). If m>0 and k + 
m>0 then g m ,h m S C 1 (]i? m i n ,oo[) with 

g' m = -mg m _ 1 , ti m = -2muh m _i. 
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Proof : The continuity assertions can be obtained using Lebesgue's domi- 
nated convergence theorem. By the same tool one can show that the func- 
tions are left differentiable with the given derivatives. Since these left deriva- 
tives are again continuous the functions are differentiable as claimed. The 
details are fairly lengthy but largely technical, and are therefore omitted. 

□ 

Using the lemma above we express p and p as functionals of U or p, 
respectively. To this end we introduce for a> —1, b> — 1 the constants 

/■i b r(a+i)r(b+i) 

c a ,b= s 1-s) ds= , 2.11 

Jo T{a + b + 2) 

where T denotes the gamma function. It will be useful later to note that 

Ca,ft-i a + b + l 

— ■ = ; , a>— 1, o>0, (2.12) 

C a ,b 

which follows from the functional relation xT(x) = r(x + l). If 4> is as in the 
lemma above and / given by the ansatz ( [2.10D then in the non-relativistic 
case we obtain 

p{r) = 2 l+3 / 2 7rc^ 1/2 r 2l g l+1/2 (U(r)), (2.13) 
p(r) = 2 l+5 / 2 7rc lA/2 r 2l g l+3/2 (U(r)). (2.14) 

For the relativistic case we substitute E 2 = (E 2 — e 2 ^) +e 2 ^ into (O) and 
use the abbreviation 

F / (e^) = e-( 2i+4 )^ i+3/2 (e' 1 ) + e -( 2 '+ 2 )^ i+1/2 (en (2.15) 

to obtain 

p(r) = 27rc l ^ 1/2 r 2l H l (e^), (2.16) 

p(r) = 2n Cl ^ 2 r 2l e^ 2l+4 ^h l+3/2 (e^). (2.17) 

A steady state of the Vlasov-Poisson system is now obtained by solving the 
equation 

C/ '( r ) = r^^J s2l+2 9l + i/2(U(s))ds, r>0, (2.18) 



which is ( |1.2| ) with (2.13) substituted in, in the case of the Vlasov-Einstein 
system one needs to solve 

e- 2A (2rA'-l) + l = le^V^r 2 ^ 2 ^), (2.19) 
e- 2X (2rp' + l)-l = 16ir 2 c W2 r 2l+2 e-( 2l+ ^h l+3/2 (e»), (2.20) 
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which is ( |1,5| ), ( |1.6| ) with (2.16) and ( |2,17j ) substituted in. It is known that 
these problems have solutions on [0,oo[, and we state the corresponding 
result for further reference; the problem of interest of course is whether 
these solutions lead to steady states with compact support and finite mass. 



Theorem 2.3 Let 6 be as in Lemma 2.k 



(a) Let Uq G K. Then there exists a unique solution U G C 1 ([0,oo[) of 
fiUf{) with U{0) = U . 

(b) Let [iq£1R. Then there exists a unique solution (A,/i) G C 1 ([0,oo[) 2 of 
fiaiQ ), ( \2JX\ ) with A(0)=0, n(fi) = i4). 



Proof : In both cases local existence on some interval [0,5] follows by a 
contraction argument, cf. [|], Thm. 3.6] for details in the non-relativistic 
case and fllC| , Thm. 3.2] for the relativistic case. In the non-relativistic 
case global existence is simple: Clearly, U is increasing. Either U < Eq on 
its maximal interval of existence, in which case U exists globally in r, or 
U{r)> Eq, r>R for some R>0, in which case p(r) = 0, r>R, and again 
U exists globally. For the relativistic case the inequality 2ra(r)<r/2 has 
to be controlled which makes the argument more involved than in the non- 
relativistic case, cf. [l(], Thm. 3.4]. □ 
Remark: In the theorem above nothing is said about the boundary condi- 
tions at infinity. However, once we know that a steady state has finite mass 
then U or [i has a finite limit at infinity, cf. Theorem |2.1| . Subtracting this 
limit from U or fi respectively and redefining Eq accordingly gives a steady 
state with the same /, but which now satisfies the boundary condition at 
infinity; the boundary condition for A follows from ( |2.7| ) if the ADM mass 
is finite. 

3 The main result 

The following theorem is the main result of the present paper: 
Theorem 3.1 Let A;,ZgIR be such that 

11 3 

k>-l, l >-^ k + l +2 >0 > k < l + 2' 

Let (j):]E m i n ,oo[— > [0,oo[ be measurable and such that<j)£Ly^ c (]E m i n ,EQ[), and 
cf>(E) = c(EQ-E) k + + 0((EQ-E) k + +5 ) asE^Eo- (3.1) 
for some Eq > E m - m , c>0, and 5>0. 
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(a) Let (f,U) be a steady state of the Vlasov-Poisson system in the 
sense that f(x,v) = (fi(E)L l with E and L as defined in (2.1 ), and 
Ue C 1 ([0,oo[) satisfies Then the steady state has compact sup- 
port and finite mass. 

(b) Let (/, A,/x) be a steady state of the Vlasov- Einstein system in the 
sense that f(x,v) = 4>{E)L l with E and L as defined in (2.k ), and 
A,/i G C 1 ([0,oo[) satisfy (\l.d[) , Then the steady state has compact 
support and finite ADM mass. 

Clearly, if <j) is as in the theorem above then it satisfies the assumptions in 
Lemma |2.2| so that this lemma and Theorem 2.3 apply. The main tool in 
the proof of Theorem 3.1 is the following lemma, which is an adaptation of 
P, Thm. 1] to our present situation: 

Lemma 3.2 Let x,y G C 1 (]0,i?[) be such that x,y>0 and 

, i s . a;+7i(r)y 
rx = a(r)y — x + x- 



l-y 2 (r)x 
ry = y[ c- (3(r)- — — 



on ]0,R[, where c>0, a, /3, 71, 72 G C(]0,R[) with oq = mf r ^ ^a(r) > 0, 
lim r ^/?/5(r) = /3oG]0,c[, 71,72 >0 ; and lim r ^ R 7i(r) = lim r ^_R72(r) = 0. Also 
let I-72 (r)x(r) > 0, re]0,R[. ThenR<oo. 

Proof : As a first step we show that there exists r* G]0,i?[ such that x(r*) > 1. 
If not, then x(r) < 1, r£]0,R[. By assumption there exists roG]0,i?[ such 
that 0(r)>O, r€[r Q ,R[. Choose K>0 such that K(l + x(r ))-y(r )>0 
and Kao > 1 + c. If for some r > ro, K(l+x(r)) — y(r) = Q then 

r(K(l + x) — y)(r) = Krx'(r) — ry'(r) 

> Kaoy(r) — Kx(r) — cy(r) 
= (Ka -c-l)y(r)+K>0. 

This implies that no such r exists, and 

K(l+x(r))-y(r)>0, re[r ,R[. 

By our assumption x < 1, and 

y(r)<K(l + x(r))<2K, re[r ,R[. 
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This implies that R = oo, and 

ry >y[c-f3{r)— — — > — — y 

for all r > n sufficiently large, cf. the assumption on (3. Integration of this 
inequality implies 



y(r)>i/(n)(r/r 1 )( <H3l °/ 2 - 



00, r — » 00, 



a contradiction. 

Thus we can assume that there exists some g]0, i?[ with x(r#) > 1. Now 
rx' > — x + x 2 = x(x — 1) (3-2) 
which implies that x(r) > 1 for all r £ [r*,i?[, and upon integration of (|3.2|) , 

Since the term in parenthesis vanishes for r = r*x{r*)/ (a?(r*) — 1) this implies 
that 

x(r*) 



'x(r*j — 1 

□ 



Proof of Theorem 3. 1 



Step 1 — The basic set-up: Consider a solution of the reduced field equa- 
tion^) as given by Theorem |2.3| . Consider the non-relativistic case first 
and define [0,R[ as the maximal interval on which U <Eq; we may assume 
that U(0) < Eq, or else the solution is trivial. If R< 00 then U(R) = Eq. If 
R = oo then Uoo = lim r _ >00 U(r) < Eq exists by the monotonicity of U. As- 
sume Uoo <Eq. Then ( |2.13| ) and the monotonicity of gi + \/2 imply that 

p(r) > 2 l+3 / 2 TTC l ^ 1/2 r 2l g l+1/2 (U oc ) = cr 21 , r > 0, 



where c>0. Then ( |l.2|) implies that U'(r) > cr 1+21 , r>0, with a different 
positive constant c, and integrating this estimate implies that £700 = 00, a 
contradiction. Thus we have the following 

Basic set-up: There exists some i?e]0,oo] such that U exists on [0,R[ with 
U < Eq on this interval, and 

lim U(r)=E . 

r— >R— 
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Analogously, we have for the relativistic case the 

Basic set-up: There exists some i?€]0,oo] such that A,/i exist on [0,R[ with 
< Eq on this interval, and 



lim e^ r) = E . 

r— >R— 



We may assume that e^ ) < Eq, or else the solution is trivial. We choose R 
maximal such that e^<Eo on [0,i?[. The non-obvious case is R = oo. By 
monotonicity, p OD = \\m. r ^n^ p{r) <Eq exists. Assume that /u tX) <£'o- Then 
again p{r) >cr 21 , r >0, with a positive constant c>0. By fl2.7|) and 



p'(r) = (l y s 2 p(s)ds) ( — J s 2 p(s)ds + 4:TTrp(r] 



>C I s 2+2l ds = cr l+2l r>Q _ 



r 2 Jo 



r 



Integration of this inequality implies that poo = oo, a contradiction. 

What we need to show in both cases is that R< 00. 
Step 2 — New variables: We introduce new variables which bring the system 
into the form stated in Lemma |3.2|. We define for the non-relativistic case 



rj(r) = E -U(r), (3.3) 

and for the relativistic case 

V (r) = lnE -fi(r); (3.4) 
recall that Eq > in the relativistic case. Now define in both cases 

x ( r ) = H*f± = JZL [ r s 2 p{s)ds, (3.5) 



rr](r) rr](r) Jo 

y{r)=^r 2 ^fl (3.6) 
p(r) 

on the interval }0,R[ with R from the previous step; note that rj,p,p>0 on 
that interval. In the non-relativistic case, rrf = —rU' = —m/r, whereas in 
the relativistic case by ( |2.9D , 

, , r]x + yp 2 /p 2 

rr] = —rp = . 

1 — 2rjx 

Thus in both cases 



, x + r )y 

rrj=-rj- —-, (3.7 

l-j 2 {r)x 
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where in the non-relativistic case 

7i = 72 = 0, 

and in the relativistic case 



7i = — o> 72 = 2r/. 
VP 2 



(3 



(3 



Now 



where 



, 47rr 2 p rrj' x + 7i(r)y 

rx = x — x — = a(r)y — x + x , 

T) l — j2{r)x 



V 



a - 



p_ 
rip' 



In the non-relativistic case we find, using ( |2.13| ), ( |2.14 ), and ( 2.12| ), that 

1 9l+3/2(U) 



a - 



1 + 3/2 r,g l+1/2 (Uy 
In the relativistic case we find, using ( p,16| ), ( |2.17| ), and ( 2.12j ), that 

1 h l+3/2 (e") 

Q. = ■ 

21 + 3 r,h l+3/2 (e") + e 2 »rjh l+1/2 (e") 
Next consider the equation for ry'. In both cases we find 

. 2rp' rp' 

ry =2yH y y. 

P P 

In the non-relativistic case Lemma U and (glf ), imply 

2/ 

r 

and 



(3. 



(3. 



= ^p - (I + l/2)2 l ^/\ Cli . 1/2 r 2l gi _ 1/2 (U) U> 
, 21 



P 



-p-pU'. 



(3. 



With the definition of r] and ( p. 7] ) this implies that 

x + 7i(r)y 



ry / = y((!W + 2)-/9(r)| 



where 



j 2 (r)xj ' 

^_^ +(2!+1) ,^|) 

P 9i+i/2{U) 



(3. 
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In the relativistic case Lemma [D| and (fn|), p7T7| ) imply 

, 21 



p =- /0 -27rq i _ 1/2 r 2/ ^(e / V 



where 



H^) = (2/ + 4)e-( 2i+4 )^ i+ 3 /2 (e^) + (4/ + 5) e -( 2i+2 )^ /+1/2 (e^) 

+ (2/ + l)e- 2 ^V 1/2 (e^), (3.14) 



and 



21 

p' = —p-(p+p)n', 
r 



which is the Tolman-Oppenheimer-Volkov equation. With the definition of 
r\ and (^) this implies that 



ry' = y((2l + 2)-(3(r)j 



where 



x + 7i(r)y 



/3 = -?7- — + 2r/ , ( 



(3.15) 



Thus both in the non-relativistic and in the relativistic case we obtain a 
system of the form which is stated in Lemma [T^, and in the next two steps 
we will show that 71,72,0;, and (3 satisfy the necessary assumptions. 
Step 3 — Application of Lemma |3. j| the non-relativistic case: In this 
case 71 = 72 = 0, cf. (|3.8|) , so these functions satisfy the assumptions in 
Lemma 3.2. To investigate the asymptotic behaviour of a and j3 we need to 
use the asymptotic expansion of <p; note that by Step 1, r/(r) — »0 for r^R. 
First of all we may assume c=l in (3.1), since this factor cancels in a and 
(3. A simple computation shows that 



/ (E -E) a (E-u) b dE = c a ^E -u) 

Ju 



a+b+1 



, u<Eq, a>— 1, b>— 1 

(3.16) 



where c a ^ is defined by fl2.11| ). Thus 

9m («) = c Km n k+m+1 + O ( v k + m + 5+1 ), v = E -u^0 + . 
Using ( $7L7\) and ( pT2| ) in ( ^T0| ) we find that 

+ 0(?7(r) 5 ), r^R-. 



(3.17) 



a(r) : 



* + !+§ 
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Since we assume U (0) < Eq, 

1 ] 1 + 3/2 (E -U(0))g l+1/2 (U(0)) > ' 

and since a > on ]0,R{ we have shown that 

inf a(r) > 
re]o,fl[ 



as required. Using ( 3.17 ) and ( |2.12| ) in ( |3.13| ) we see that 

/ 3(r) = -(fc + ; + 5/2) + (2/ + l) fc ||+ / 3 2 /2 +0( ?? (r) 5 )^fc + ; + ^ r^R-, 

and by our assumptions on k and I this limit lies in the interval ]0,2Z + 2[ 
as required. Applying Lemma 3J2 we find that R<oo, and the proof of the 
theorem is complete in the non-relativistic case. 

Step 4 — Application of Lemma \3. 4 , the relativistic case: By (|3.9| ) and since 
i](r)^0 as r^R—, 72 is as required. We use fl3,l|) to compute the asymp- 
totic behaviour of the various other quantities as r— > R— . To this end we 
first observe that ( |3.1|) implies that 



<j ) {E) = dE{E 2 -E 2 ) k + +EO({E 2 -E 2 ) k + +5 ) as E^E -, 

which is more suitable for computing the integral h m . Again we assume 
without loss of generality that d = 1. Using the relation 

/ E{E 2 - E 2 ) a (E 2 - u 2 ) b alE = -c a , b (E 2 - u 2 ) a+b+ \ u < E , a > -1, b > -1, 

Ju 2 

we find that 

M^) = \c k , m (E 2 -e 2 n k+m+1 + 0((E 2 -e 2 n k+m+S+1 ) 

= \c Km e k+m+1 + 0{e k+m+5+1 ), e = E 2 -e 2>1 ^0 + . (3.18) 
In the following we need to observe that 

' e + 0(e 2 ), e^0+ (3.19) 



■•~2E 2 
and 



E 2 e- 2 ^ = l + 0(e(r)), r^R-. (3.20) 
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Using ( gUp , ( |3l9| ), ( pO|) , and pl2| ) in ( ^Tl| ) we find that 



a(r) 



1 



2 -2ju(r) C k,l+3/2 



+ 0(e(r) 



1 



+ 0(e(r) 5 ), r^i?-. 



2Z + 3 u " c M+1/2 ' v ~ v ' ' k + l + h/2 

Since a > on ]0,R[ and lim r ^o+o ; ( r ') > we find that 

inf a(r) > 
re]o,R[ 

as required. Since by fl3,9|) 71 = ap/p, ( 2.16 ) and ( 2,17 ) together with Q3,18|) 
imply that 71 (r)-»0 as r-> R-. 

It remains to examine the function [3. Using ( 3.18 ) in ( |2.15| ) and ( |3.14| ) 
we find 

H,(ef) = I Cfc ^ 1/2 (2/ + l) e - 2 ^ e ^+ 1 /2 + 0(e ^ + i/2 )) 
and using this together with (|3.19| ), ( |3.20 ), and (|2.12| ) in ( 3.15[) we obtain 
m = 0(e(r)) _ a -i (r) + 2(2Z + 1) 6 Cfc ''- 1/2 + 0(e*(r)) 

^0 c fe,/+l/2 

= -(£; + / + 5/2) + 2(£; + / + 3/2) + 0(e<V)) 

-» A: + / + ^, r^R-, 

and by our assumptions on & and I this limit lies in the required interval 
]0,2Z + 2[. Applying Lemma 3^ completes the proof in the relativistic case. 

□ 



Remark: We note that the application of Lemma |3.2| provides an explicit 
upper bound on the radius R of the spatial support in terms of a point r* > 
at which 

m(r*) 



namely 



x(r*) ■ 



R<r* 



>1 



and this upper bound 



x[r*) — 1 

Such a point must exist by the proof of Lemma | 
may be useful in numerical work on steady states. 

We have shown that the spatial support of a steady state as in Theo- 
rem 3.1 is compact, but a bound on v over the support of / follows by the 
boundedness of E from above together with ( p.ip or Q2.2] ) and the bound- 
edness of U or fi respectively. 
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4 Examples and final remarks 



The following is a list of types of static spherically symmetric models built 
of self-gravitating collisionless matter, non-relativistic and relativistic, to be 
found in the literature. In this list inessential multiplicative constants are 
omitted. 

(NR1) Polytropic solutions of Vlasov-Poisson system. The distribution func- 
tion is of the form (Eo — E) k L l for E<Eq and zero otherwise. For 
appropriate k and I — cf. the introduction — the existence of solutions 
with finite radius was proved in Theorem [O] applies, but only to 
a subclass of these, and we comment on this fact below. 

(NR2) King models, cf. [||, p. 232]. The distribution function is of the form 
e E °~ E — 1 for E<Eq and zero otherwise. Since 

e E -E _ 1 = {Eq _ E) + 0(( £ o _ £ )2)) E ^ £ q _ 



Theorem 3.1 applies to these models. 

(NR3) Woolley-Dickens models, cf. [§, p. 235]. The distribution function is 
of the form e E °~ E for E < Eq and zero otherwise. Since 

Theorem fO applies. 



(NR4) Wilson models, cf. |2|, p. 235]. The distribution function is of the form 
e E °~ E — 1 — (Eq — E) for E < Eq and zero otherwise. Since 

e E -E_ 1 _ {Eo _ E) = ^ Eo _ E) 2 + 0{{Eo _ E) 3 ) ^ £^£ _ 

this model is outside the range of the present approach. 



(Rl) Truncated Maxwell-Boltzmann models, cf. [ pi| p. 59]. The distribution 
function is given by an expression formally identical to that of the 
Woolley-Dickens models (of course E has a different definition in the 
two cases), and Theorem [O] applies. 

(R2) Power-law models, cf. |2l|, p. 68]. The distribution function is given 
by [(E/Eq) 2 ]- 5 [1-(E/Eq) 2 } 5 for E<E and zero otherwise. Since 

[(E/Eq) 2 ]- s [1-(E/E ) 2 ] s = ^(E -E) s + O((Eq-E) 2S ), E^Eq- 



Theorem 3.1 applies if -1/2 < 5 < 3/2. 
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(R3) Polytropic models, cf. [21, p. 68], The idea here is the follow- 
ing. Polytropic solutions of the Vlasov-Poisson system in the sense 
of (NR1) with ^ = correspond to self- gravitating fluid models with 
p = p(n+i)/n^ w j iere p i s ^he density and p is the pressure. This is due 
to the relation between kinetic and fluid models mentioned briefly in 
the introduction and discussed in more detail below. In general rela- 
tivity it is also possible to look for a collisionless model with a given 
equation of state but there is no unique generalization of the poly- 
tropic case. In [^] Tooper considered different possibilities and their 
relationships. One possibility is to consider the equation of state which 
is formally identical to the non-relativistic polytropic one with p being 
interpreted as the energy density. If, on the other hand, p is interpreted 
as the mass density a one-parameter family of relations between pres- 
sure and energy density is obtained. Some authors, e.g. Packer ell |||], 
have considered the problem of producing corresponding distribution 
functions. Numerical calculations indicate that this is sometimes pos- 
sible but not always. We can prove the relative statement that when it 
is possible (and the polytropic index n is restricted appropriately) the 
radius of the configuration is finite. Note that the in the case where 
depends on E alone the asymptotic behaviour of goes into the proof 
of Theorem 3.1 only via the asymptotic behaviour of the equation of 
state of the corresponding fluid model. In the first case the equation 
of state is directly in the correct form. In the second case the rela- 
tion between pressure and energy density is p = Cp n l <yn+l ' > + np for a 
positive constant C. Then for small p 

p = c-( n + 1 )/ n /9 ( n + 1 )/ n (i + p l l n ) 

Next we want to state two observations which may help to understand 
the relevance of the upper bound k < 1 + 3/2 that we required in our ansatz. 
We restrict ourselves to the non-relativistic case. The first observation is 
based on the correspondence of steady states of the Vlasov-Poisson system 
with steady states of the Euler-Poisson system. If (f,U) is a steady state 



of the Vlasov-Poisson system with f(x,v) = <j)(E) with (f> of the form (3.1), 
then (p,p,U) satisfy the Euler-Poisson system: Since the solution is static, 
i. e., 

>3 



u(x) = [ vf(x,v)dv = 0, xGJR 3 
p[x)J 



p{x) 

the equation p' = —pU' is all that remains of the Euler equations, and this 
is ( |3.12 ) in the isotropic case 1 = 0. Since the functions g m in Lemma [T^ 
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are strictly decreasing on their support, p can be written as a function of p, 
i. e., one obtains the equation of state 



V = ^1/2(02(5-1/2) 1 (pj) =C3p n ™ 1 +0(p n « 1+e ), p^0+ 

where n = fc + 3/2, ci,C2,C3 are positive constants which depend on the con- 
stants appearing in ( 2.13| ) and ( 2.14 ), and e>0 depends on 5. The critical 



value k = 3/2 corresponds to n = 3, and we now give an example of a steady 
state of the Euler-Poisson system, which has an equation of state of poly- 
tropic form with 3 < n < 5 for p small, and with unbounded support and 
infinite mass. It is based on a well-known explicit singular solution (cf. |J, 
p. 89). 

Example: Let n be a real number with 3 < n < 5 and let M be the set of 
functions p{r) from [0,oo[ to itself satisfying the following conditions: 

1. p is a smooth (C°°) function of x = r 2 

2. p'(r)<0 for r>0 

3. p"(0)<0 

2n 

4. p(r) = r n-i forr>l 

This set of functions is convex. We claim that there is a function p in A4 
with J l r 2 p(r)dr = ^E^- To see this, note first that the function ps( r ) = 

2n 

r n-i satisfies this condition. However, it does not belong to M. Let 
P- be an element of M which is everywhere less than or equal to ps- It 
is clear that such a function exists. It satisfies the condition that I_ = 
Jq 1 r 2 p^{r)dr <^E^- It is also clear that there exist functions p + £M such 
that I + = Jq r 2 p + {r)dr is as large as desired. In particular p + can be chosen 
so that I + > ^Ejj- There exists a A in the interval [0, 1] such that AJ_ + (1 — 
A)i+ = ^5^. Then p = Ap_ + (1 — \)p+ is the desired function. 

Define m(r)=47r Jq s 2 p(s)ds, where p is the function just constructed. 
Let ms be the function constructed in the corresponding way from ps- By 
construction Jq 1 s 2 p(s)ds = Jq 1 s 2 ps(s)ds. Thus m(l) =ms(l). It follows that 
m(r) = ms(r) for all r > 1. Now define p(r) = s~ 2 m(s)p(s)ds and ps anal- 
ogously. Then for r > 1 the functions p and ps are equal. Since both p and 
p are strictly monotone for r > it is possible to write p = f{p) for a smooth 
function / for p in the interval (0,p(0)). Moreover, df /dp is strictly positive. 
Both p and p are smooth functions of x for x > 0. If it could be shown that 
dp/dx and dp/dx are non- vanishing at x = then it would follow that / has a 
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smooth extension to the interval (0,oo) with everywhere positive derivative. 
This is equivalent to showing that d 2 p/dr 2 and d 2 p/dr 2 are non-vanishing 
at r = 0. The first of these two quantities is non- vanishing by assumption. 
The other can easy be computed to be equal to — ^-p 2 (0). The conclusion is 
that the functions p and p are related by an equation of state satisfying all 
the usual conditions, i.e. that for p>0 the conditions p>0 and dp/dp>0 
are satisfied everywhere. 

For r > 1 the function m(r) can be calculated explicitly, as it is equal 

to ms(r). The result is m(r) =4ir ^3 r n ~ 1 . In the same way p(r) can be 

( n _]\2 -2(n + l) 

calculated explicitly for r>l with the result p(r) = 4ir 2 ( n " i)(n-3) r n_1 
Comparing with the expression for p shows that the equation of state is of 
polytropic form with index n for p< p(l). 

The example is not dependent on the condition n < 5. That restriction 
was only made because that is the interesting case where it is known that 
there are some solutions with finite mass and radius, namely those with 
polytropic equation of state. Note also that it was shown in [18| that there 
exist some solutions of the Euler-Poisson and Euler-Einstein equations which 
have an equation of state which is asymptotically like a polytropic one with 
3 < n < 5 (but not exactly polytropic) which have finite radius. This was 
done by a perturbation argument with the attendant disadvantages. The 
above example serves to show that the strategy applied by Makino^] cannot 
be directly extended to the case n > 3 and strongly suggests that the methods 
of this paper cannot be modified to cover cases with k > 3/2. The polytropic 
steady states are in some sense structurally unstable. 

The second observation is the following: Consider for k > the energy- 
Casimir functional 

defined on the set 

F M = {f£L 1 nL 1+1/k (M 6 )\f>0, J J fdvdx = M, J Jv 2 fdvdx < 00}, 

where Uf denotes the potential induced by / and M > is a prescribed 
constant. One can show that V is bounded from below on Tm if k<3/2, 
whereas this functional is not bounded from below on this set if k> 3/2. In 
the first case one can then show that the functional actually has a minimizer, 
at least if one restricts the set Tm to spherically symmetric functions, and 
this minimizer is a steady state of the Vlasov-Poisson system of polytropic 
form which is dynamically stable in a well-defined sense, cf. [||, |6| . Again we 
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see that there seems to occur a loss of stability if one crosses the threshold 
fc = 3/2, i. e., n = 3. In the Casimir part of the functional T> the expression 
ji+i/k can re pi acec | D y m ore general functions Q(f), provided Q satisfies 
growth conditions both for small and for large values of /, which have to 
satisfy the same restrictions as k above. This leads to stable steady states 
not necessarily of polytropic form. 
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